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Abstract 

We study the electron transport in open quantum-dot systems described by the interacting 
resonant-level models with Coulomb interactions. We consider the situation in which the quantum 
dot is connected to the left and right leads asymmetrically. We exactly construct many-electron 
scattering eigenstates for the two-lead system, where two-body bound states appear as a con¬ 
sequence of one-body resonances and the Coulomb interactions. By using an extension of the 
Landauer formula, we calculate the average electric current for the system under bias voltages in 
the first order of the interaction parameters. Through a renormalization-group technique, we arrive 
at the universal electric current, where we observe the suppression of the electric current for large 
bias voltages, i.e., negative differential conductance. We find that the suppressed electric current 
is restored by the asymmetry of the system parameters. 
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I. INTRODUCTION 


In the last two decades, much progress has been made in the experimental studies of 
the electron transport in nanoscale devices^ - ! In the systems smaller than the coherent 
length, quantum effects are observed in the electron transports^. In order to analyze it 
beyond the linear-response regime theoretically, we need to treat nonequilibrium steady 
states realized in open quantum systems. The Landauer formula/’ 8 enables the calculation of 
the electric current flowing across nanoscale samples in noninteracting cases, which indicates 
that the nonequilibrium steady states are scattering states in the open quantum systems. 
Indeed, the transport properties such as the electrical conductance and the electric-current 
noise are determined by the scattering matrix- - — . To investigate interacting cases, the 
Keldysh formalism of the nonequilibrium Green’s function has been developed- 12- - 6 . It has 
provided a standard tool for the study of the Kondo effect measured as a conductance peak in 
semiconductor quantum dots (QDs)^ t18 . On the other hand, we have proposed an extension 
of the Landauer formula to interacting cases^ 1 ^, and have shown that the scattering states 
are essential in the interacting cases as well. 

The interacting resonant-level model (IRLM) is one of the standard testbeds for such stud¬ 
ies of the open QD systems with interactions. The original IRLM, which consists of a single 
impurity coupled to a conduction band, was introduced for studying the Kondo problem in 
equilibrium systems^. Recently, the IRLM with two external leads has been employed as a 
minimal model of open QD systems with Coulomb interactions; it plays an important role in 
verifying the theoretical approaches such as the nonequilibrium Bethe-ansatz approach- 22 , the 
perturbation theory with the numerical renormalization group^, a new method called impu¬ 
rity conditions^ and the time-dependent density-matrix renormalization-group method^. 

A remarkable feature of the two-lead IRLM is the appearance of negative differential 
conductance , that is, the suppression of the electric current due to the Coulomb interaction 
for large bias voltages. Clearly, this is a phenomenon out of the linear-response regime. To 
see the feature and to compare the results obtained by different approaches, the universal 
electric current characterized by a single scaling parameter Tk is useful- 24 ’^S^. Indeed, it is 
found that, for large bias voltages V, the universal electric current shows a power-law decay 
(I) oc (V/Tk)~ u ^ with the parameter U of the Coulomb interaction^ 4 ’—. 

In the previous papers^^, we proposed an extension of the Landauer formula with 
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many-electron scattering eigenstates. We considered the two-lead IRLM with linearized 
dispersion relations and gave exact many-electron scattering eigenstates in explicit forms. 
This is in contrast to the previous studies^ - — of the scattering problems for other QD 
systems, which include integrals or matrix inversions. The explicit iV-electron scattering 
eigenstates enabled us to calculate the quantum-mechanical expectation value of the electric 
current, which we called the iV-electron current. By taking the electron-reservoir limit 
N —> oo of the IV-electron current, we obtained the average electric current for the system 
under finite bias voltages. It is clear that the way of realizing the nonequilibrium steady 
states in our extension of the Landauer formula is different from the Keldysh formalism^ - —. 
By employing a renormalization-group technique with the Callan-Symanzik equation 24,2 -, 
we arrived at the universal electric current in the first order of the interaction parameter 
U. We found that the negative differential conductance of the universal electric current is 
characterized by the same scaling parameter Tk as that obtained by other approaches^^L^k 


We remark that the apparent inconsistency in Ref. 19| pointed out in Ref. 
the level of the universal electric current 43 . 


32 


is removed in 


In the present paper, we study the two-lead IRLM in which the QD is connected to 
the two external leads asymmetrically. The effect of the asymmetry of the QD systems 
is observed in experiments. For example, in semiconductor QDs^, the asymmetry of the 
lead-dot couplings causes breaking of the unitary limit of the Kondo effect, which is theoret¬ 
ically understood in the linear-response regime^ 1 ^. In the present study, we investigate the 
effect of the asymmetry on electron transport out of the linear-response regime. One of the 
theoretical difficulties of the asymmetric cases is that the even-odd transformation, which 
maps the two-lead IRLM to two single-lead systems^>22,> 22 ) does not work. The application 
of the Bethe-ansatz approach- 22 has been restricted to the cases in which the even-odd trans¬ 
formation works. The construction of exact many-electron scattering eigenstates for such 
pure two-lead systems is established for the first time in this paper. 

Through the extension of the Landauer formula and a renormalization-group technique, 
we obtain the universal electric current for the asymmetric cases in the first order of the 
Coulomb-interaction parameters U\ and U 2 - The universal electric current is characterized 
by the two renormalized parameters Ti and T 2 . The sum Tk = T x + T 2 provides a scaling 
parameter for the bias voltage V, which is similar to the symmetric cases. Our universal 
electric current has the same functional form as that obtained by the renormalization-group 
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approach^A^, although there is a difference in our calculation of the renormalized parameters 
T\ and T 2 . As we will point at the end of Section HV1 this leads to a critical difference in the 
predicted behavior of the universal electric current. The suppressed electric current due to 
the Coulomb interaction is restored by the asymmetry of the system parameters. To clarify 
the relation between the asymmetry of the system parameters and the restoration of the sup¬ 
pressed electric current, we introduce the asymmetry parameter 5 = {U\ — U?)(Ti — T 2 )/(2Tk) 
taking the value in the range 0 < S < U with the average interaction U = (U\ + T/ 2 )/2. In 
fact, in the first order of U\ and C/ 2 , the power-law decay of the universal electric current is 
given by (I) oc (V/T^)~^ u ~ 5 ^ w , which indicates that the restoration of the suppressed elec¬ 
tric current occurs with both asymmetric Coulomb interactions and asymmetric lead-dot 
couplings. The restoration was reported to happen even for symmetric lead-dot couplings 
in Refs. 
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and 
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but we presume that this was due to higher orders of the interaction Ug 


in the renormalized parameters T\ and T 2 . 

The exact many-electron scattering eigenstates tell us much about the transport proper¬ 
ties of interacting electrons in the open QD systems. We notice that the scattering processes 
in which the set of wave numbers of incident plane waves is not conserved are essential in 
interacting cases. The explicit form of the scattering eigenstates indicates that, due to the 
Coulomb interactions, the incident plane-wave states are partially scattered to many-body 
bound states that decay exponentially as the electrons separate from each other. Indeed, the 
two-body bound states appear in the two-lead IRLM; each term of the //-electron scattering 
states is characterized by the configuration of the two-body bound states. We can under¬ 
stand the origin of the negative differential conductance in the two-lead IRLM in terms of 
the formation of two-body bound states. Such many-body bound states are also found in 
other open QD system s 34 ’ 35 . 

The present paper is organized as follows. In Sec. Ill Al we introduce the two-lead IRLM 
with asymmetric interactions. In Sec. IIIBI the extension of the Landauer formula^ is 
described in a general setting. In Sec. IIII1 we present the construction of the exact one- 
and two-electron scattering eigenstates whose incident states are free-electronic plane waves 
in the leads. We also give the //-electron scattering eigenstates in the first order of the 
Coulomb-interaction parameters. In Sec. lIVl through the extension of the Landauer formula, 
we calculate the average electric current for the system under finite bias voltages. We 
obtain the universal electric current by dealing with the divergences in the average electric 


4 








current with the Callan-Symanzik equation. As a result, we observe the negative differential 
conductance and the restoration of the suppressed electric current. Section [V] is devoted to 
concluding remarks. 

II. MODELS AND FORMULATION 

A. Interacting resonant-level models 

We consider the open QD system described by the IRLM of spinless electrons. It consists 
of a QD with a single energy level and two external leads of noninteracting electrons. The 
arrangement of the QD and the two leads is illustrated in Fig. |T| We assume the situation 
in which the QD is connected to the two leads asymmetrically. 

The Hamiltonian is given by 

n Id 

H = ^ / dx c\(x)--^—cg(x) + ^ (tgcl(0)d + t* e d^cg(0)) + e d d)d 

/ L 1 CIOC 

t = 1 , 2 j ~2 1 = 1,2 

+ U t c\(0)c e (0)<fld. (1) 

€= 1,2 

Here cj(x) and cg(x) are the creation- and the annihilation-operators of an electron in the 
lead £, while d I and d are those on the QD. The first term corresponds to the kinetic energy 
of electrons in the leads, where L stands for the length of the two leads to be eventually 
taken to infinity. The second term expresses the tunneling between the leads and the QD, 
where the parameter tg(E. C) is the transfer integral. We assume a single energy level e d (G M) 
on the QD, which corresponds to the third term. The fourth term describes the Coulomb 
interaction between the two electrons at the origin x = 0 of the lead i and on the QD. The 
parameter Ug{> 0) is the strength of the Coulomb repulsion. 

We focus on the electrons with positive velocities in the vicinity of the Fermi energy 6p of 
each lead and linearize the local dispersion relation to be e{k) = Vp{k — kp ) + ep under the 
assumption that the other parameters \tg\, |e d | and Ug are small compared with the Fermi 
energy ep^i. For simplicity, we have set vp = 1, kp — 0 and ep = 0 in Eq. (|T[) . Then, as is 
indicated in Fig. [TJ an electron coming from the part x < 0 of the lead 1 is scattered at the 
QD to the parts x > 0 of the two leads. 

In constructing the scattering eigenstates, we treat the system as an open system in 
the limit L —> oo of the two leads. In the sprit of the original Landauer formula 7 ^, we 


5 



FIG. 1: The two-lead interacting resonant-level model. Electrons flow only upward due to the 
linearized dispersion relations of the leads. 

suppose that the infinite two leads can substitute for electron reservoirs that are in the Fermi 
degenerate states of noninteracting electrons. We assume that the electrons emitted from 
the electron reservoir into the lead i follow the Fermi distribution function / Wi( a e (E) = 1/(1 + 
e +(B-/+) w ith the chemical potential /i^ and the inverse temperature fa. We are interested 
in the nonequilibrium steady state realized between the large two electron reservoirs in the 
cases /ii ^ /j ,2 of different chemical potentials. 

We adopt a standard definition of the electric-current operator as 

/ = ah - (1 - a)I 2 , l£ = i(tec\(0)d-t* e dfae(0)). (£=1,2) (2) 

For an arbitrary eigenstate \ij)) of the Hamiltonian H , the expectation value fa>\I\^) does 
not depend on the parameter a since the relation {fa\I]\^) = — (^|/ 2 | , 0) holds. In what 
follows, we choose the parameter a = \t 2 \ 2 /t 2 with t = \/|fi| 2 + |^ 2 1 2 for the convenience of 
calculations. 


B. Extension of the Landauer formula 


Our purpose is to study the average electric current flowing across the QD of the two-lead 
IRLM beyond the linear-response regime. The extension of the Landauer formula, which 
was proposed in Refs. 


19 and 


2d. consists of the following three steps: 


(i) Construction of many-electron scattering eigenstates whose incident states are free- 
electronic plane waves in the leads; 

(ii) Calculation of the quantum-mechanical expectation value of the electric current with 
the many-electron scattering eigenstates; 
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(iii) Calculation of the statistical-mechanical average of the electric current by assuming 
the equilibration of electrons in each electron reservoir. 

The many-electron scattering eigenstates constructed in the step (i) are characterized by the 
wave numbers of the incident plane waves. We note that they are essentially different from 
the Bethe-ansatz eigenstates^ l22 . ; whose incident states are not free electronic but include 


the effect of interactions. The Bethe-ansatz result in Ref. 


22 did not agree with results of the 


previous works^^^ while our results agree with them. For the wave numbers {ki ,..., k^} 
of the IVi-electron incident plane wave coming in through the lead 1 and {hi,..., h^} of the 
W-electron incident plane wave coming in through the lead 2, we express the IV-electron 
scattering eigenstates by | k,h) = |fci,..., k^, hi ,..., h^ 2 ) with N — Ni + N 2 . In the 
step (ii), we calculate the expectation value (k,h\I\k,h)/{k,h\k,h) of the electric-current 
operator /, which we call the IV-electron current. This calculation is practically carried 
out by using the explicit IV-electron scattering eigenstates. In the step (iii), we take the 
limit Ng, L —> oo of the IV-electron current by assuming that the wave numbers fc* and hi 
of incident plane waves follow the Fermi distribution of each electron reservoir. We call 
the limit an electron-reservoir limit. Clearly, the reservoir limit corresponds to taking the 
statistical-mechanical average of the electric current for all the incident states that follow the 
Fermi distributions. In general, the electrons scattered at the QD are in many-body states 
including the effect of interactions. We assume that such many-body states are completely 
equilibrated to the Fermi degenerate state of free electrons in each election reservoir before 
being re-emitted towards the QD, which is the main assumption of the extension of the 
Landauer formula. We shall see for the two-lead IRLM that, since the IV-dependence of the 
IV-electron current appears only in the upper bounds of the sums on wave numbers {ki} 
and {hi}, we can take the reservoir limit by replacing the sums with the integrals on k and 
h with the Fermi-distribution functions / Mlil a(/c) and / M 2 ,^(h). 

The way of realizing the nonequilibrium steady states in the extension of the Landauer 
formula is different from that in the Keldysh formalism-^ ——. I n Q ur extension of the Lan¬ 
dauer formula, we first construct the IV-electron scattering eigenstates for finite N without 
the information of the equilibrium states in the electron reservoirs. After the calculation of 
the IV-electron current, we take the reservoir limit to consider the nonequilibrium steady 
state. In the Keldysh formalism, on the other hand, the Green’s functions or the density 
operator describing the nonequilibrium steady states are obtained by adiabatically turning 
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on the perturbative terms for the initial nonperturbative steady states of infinite number of 
electrons. 

Our approach is also independent of Hershheld’s bias-operator approach, which constructs 
the density operator of the nonequilibrium steady states directly from one-electron held op¬ 
erators in the framework of the quantum held theory^ - — . We remark that the construction 
of the density operator through the bias-operator approach has not been established ana¬ 
lytically in interacting cases^L^S except for the Toulouse limit of the Kondo mode l 39-40 . 


III. MANY-ELECTRON SCATTERING EIGENSTATES 

A. One-electron cases 

The linearization of the local dispersion relations of the leads enables us to construct 
exact scattering eigenstates. First, we consider the one-electron cases. The one-electron 
scattering eigenstates are given in the form 

\E) = ( [ dx Y dmixylix) + ed^ | 0 ), (3) 

771 = 1,2 

where |0) is the vacuum state satisfying q(x)| 0) = d|0) = 0. The eigenfunctions g m (x) and 
e are determined by the coupled Schrodinger equations: 

CiiL -E) 9 ™^ + t meS(x) = 0, [m — 1 , 2) 

(e d - E)e+ Cfcf0) = 0. (4) 

771 = 1,2 

It is readily found that the eigenfunction g m (x) is discontinuous at x = 0 and the matching 
condition at x = 0 is obtained by integrating the first equation in Eqs. ([4]) around x = 0 as 

g m ( 0+) - .^(O—) + 1 t m e = 0. (5) 

Since the value £ 777 ( 0 ) is not determined by the Schrodinger equations, we assume £ 777 ( 0 ) = 
(£ 777 ( 0 +) +£ 777 ( 0—))/2 from physical intuition. 

To employ the Landauer formula, we need the scattering eigenstates whose incident states 
are a plane wave in the lead 1 or the lead 2. For the incident plane wave with the wave number 
k in the lead £, we consider the solution g^ k {x) and e£* of the Schrodinger equations (0} 


satisfying 


1 


9 %( x ) = for x < °> 


( 6 ) 


where S m e is the Kronecker delta. We refer to Eq. fl 6 ]) as scattering boundary conditions. 
The solution with energy eigenvalue E = k is given by 


9m,k( x ) = y=(<W ~ it m V2ne^0(x))e lkx , 

e W = _i_5_ (7) 

v^fc-e d + ir’ [) 

where 6{x) is the step function and T = (|ti| 2 + |^2 1 :: )/2 is the level width of the QD. By 
inserting them into Eq. (j3J), we obtain the scattering eigenstate | k\£) whose incident state 
is a plane wave with wave number k in the lead £. 

The one-electron scattering eigenstates \k,£) are normalized on the (5-function as 
(fc; £\k'; l') = Se£>5(k — k') in the limit L —> oo. In the calculation of quantum-mechanical 
expectation values of physical quantities with the scattering eigenstates | k]£), we need to 
restore the length L of the leads in order to regularize the square norm as (/c; £\k; £) = L/(2 tt). 


B. Two-electron cases 


We next consider the two-electron cases as the simplest example of the interacting cases. 
The form of the two-electron scattering eigenstates is given by 



Here we impose the antisymmetric relation gi m (x 1 ,^ 2 ) = —gmi{x 2 ,aq). The eigenvalue 
problem H\E) = E\E) leads to the coupled Schrodinger equations: 


1 \dxi 
1 d 


d d 
+ 


dx<) 


E)gi m {x i,x 2 ) + t m ei(x 1 )5(x 2 ) -ti5(x i)e m (x 2 ) = 0, 


+ e d + Ui5(x) - E^Jei(x) + ^ = 0. 


m= 1,2 


(9a) 

(9b) 


In the previous works^^ for the symmetric case U\ = f/ 2 , we employed the even-odd 
transformation that maps the two-lead IRLM to two single-lead systems. However, since 
the transformation does not work for the asymmetric cases U\ 7 ^ U 2 , we deal with the 
two-lead IRLM directly. 
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We present a construction of the exact two-electron scattering eigenstates, which is an 
extension of the previous one^>2Mi First, we derive three important relations from the 
Schrodinger equations (15a|) and f[9bjl . The eigenfunction gi m (x i,x 2 ) is discontinuous at Xi = 
0 and x 2 = 0, while ei(x ) is discontinuous at x — 0. The matching conditions at the 
discontinuous points are given by 


gim{x , 0+) - gim(x, 0—) + i t m ei(x) = 0, (10a) 

e z (0+) - e*(0-) + iUiei(0) = 0, (10b) 


which are obtained by integrating the Schrodinger equations (l9a|) and (l9b]l around the dis¬ 
continuous points. Since the values of the eigenfunctions at the discontinuous points are not 
determined by the Schrodinger equations, we assume 


gim(x, 0) = ^ (gim(x, 0+) + gi m (x , 0-)), (11a) 

ei( 0) = ^(e t ( 0+) + ej(0-)) (lib) 


in a way similar to the one-electron cases. By applying Eqs. (1 10 aft and (11 lap to Eq. (l9b|) for 
x ^ 0, we have 


1 d 
i dx 


fd - iT - E\ei(x) = -FCgimfeO-)- 


( 12 ) 


Given functions gi m (x, 0—), (m = 1, 2), we obtain the general solution for ei(x) as 

ei(x) = Ci^ E ~^ x - i^C [ X dze< E - e « +lT ^g lrn (z, 0-), 

_ J xa 


(13) 


where Ci is the integration constant and xq is chosen as xo = — oo if x < 0 and xo = 0 
otherwise. On the other hand, by applying Eq. (Ill bll to Eq. (llObp . we have the matching 
condition 


(l + ^)e,(0+) = (l - jtl,)e,(0-). (14) 

The equations (llOajl . (TT3]l and (IT4|) are the relations that we need. 

Next, we demonstrate how to construct the two-electron scattering eigenstates by the 
repeated use of the three equations (llOal) . (TT3h and (fT4lh Let us consider the situation in 
which one electron with wave number k\ is coming in through the lead l\ and another with 
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(*i,40 

All,(12) 

All,(21) 

Al2,(12) 

Al2,(21) 

A 2 2,(12) 

A 22 ,( 2 i) 

(1,1) 

1 

-1 

0 

0 

0 

0 

(1,2) 

0 

0 

1 

0 

0 

0 

(2,2) 

0 

0 

0 

0 

1 

-1 


TABLE I: The coefficients Ai m)P of the incident plane-wave states for (£i,£ 2 ) = (1,1), (1,2) and 
( 2 , 2 ). 


k ‘2 is coming in through the lead £ 2 - We construct the eigenfunctions with energy eigenvalue 
E — k\ + k ‘2 that satisfy the scattering boundary conditions 


gim(x 1 , X 2 ) = ^Y Am, P e l{kp i X1+kp > X2) (15) 

P 

for xi,x 2 < 0. Here P = (Pi,P 2 ) is a permutation of (1,2) and the coefficients Ai m ^ P 
are given by Ai m: p = sgn(P)5^ p <5^^ with the signature sgn(P) of the permutation P. 
The coefficients Ai m)P are explicitly listed on Table HI Beginning from the incident state 
ge„i(xi , X 2 ) in the region x\ < x 2 < 0, we connect it to the other regions through Eqs. fllOajl . 
(TT3]1 and (TT4|) . By inserting gi m (x, 0—) into Eq. (IT3j) . we have 


ei(x) 



V A lm ,p 


P,m 


(16) 


for x < 0. Here we have set xo = —00 and have taken Ci — 0 to avoid the divergence as 
x —$■ — 00 . The function gi m (x, 0—) is connected to gi m (x, 0+) by Eq. (llOaft with Eq. (TT6h as 

9 lm(x, 0 + ) = ^ Y A ln,pe lkPlX ($mn ~ lt m V^7Te^), (17) 


which leads to 


7 ^2) 



J2 A ‘».p e ' kF ' n sA Pa (*) 


(18) 


for xi < 0 < x 2 . Recall that g„‘\{x) is the one-electron scattering eigenfunction in Eqs. (JTJ) . 
Again, by inserting gi m (x, 0—) into Eq. (fT3lh we have 


d(x) = C[A E - eA+lV)x 


Y A ™n,pgl% 2 ( x ) e( k£ 

P,m,n 


(19) 
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for x > 0. Here we keep the first term with the integration constant C[ since the term is 
not divergent as x —> oo. By inserting gi m (0—,x) and e/(x), ( x > 0) into Eq. (110 aft with the 
antisymmetric relation gi m (xi,x 2 ) = —g m i{x 2 ,x i), we have 


gim(x i,X 2 ) = XI A m,P 9 ^ Pl {Xl)g^ kp 2 {x 2 ) + ltjC^e 1 ^ lT "> x ^ E ^+1^2) ( 20 ) 

P,r,n 

for 0 < Xi < x 2 . Finally, by inserting Eqs. (TTH and (Til?]) into Eq. (|T1|) . we determine the 
integration constant C[ as 

C, =-}= m (S lh e^ - 6 u ,e <f) (21) 

with ui — 2f/;/(2 + if/;). Thus the two-electron scattering eigenfunctions satisfying the 
scattering boundary conditions (TT5|) are obtained as follows: 

9imMk 2 ( x 1^2) = 9iil( x i)9m%( x 2) ~ 9iM( x i)9m%(x2) 

+ i2)e lEx2 6(x2i)d( Xl ) - t m u l Z^l(x 2 i)e lExi e(x 12 )e(x 2 ), 

e SS ( x ) = HI ( x ) e fc 2 2) - HI ~ (-^^(x), (22) 


where x t] = x^ — x 3 and 


7 (^ 2 ) ( T \ - 

Zy m , fcifc 2 V X / — 


y/2ir 




(23) 


Here, on the left-hand sides of Eqs. (I22|) and (|23|) . we write the wave numbers k\ and k 2 and 
the superscripts i\ and l 2 of the leads explicitly. 

Each term of the eigenfunctions in Eqs. (]22l) is interpreted as follows. The first two terms 
correspond to the two-electron scattering eigenfunctions of the noninteracting cases, which 
are given by the Slater determinant of the one-electron scattering eigenfunctions in Eq. flZj). 
The effects of the interactions appear in the terms with the function (x). They are 

interpreted as two-body bound states since they decay exponentially as the two electrons 
separate from each other. For example, the third and the fourth terms in the eigenfunction 
9imMk 2 ( Xl ’ x2 ) Eq. $22]) are rewritten as 

t i u mZ ( J i 'ti 2 (xn)e' El2 0(x 2 i)0(x I ) - t m u,Zl‘^l(x 21 )e' Ex '9(x u )0(x 2 ) 

= (O)0(z 21 ) - 0)9(i 12 )) e Wf-'^- r )'"-“'+' B£1 ? 2 e(x 1 )e(x 2 ). (24) 


The binding length of the two-body bound states is given by 1/T where T is the level width 
of the QD. It should be emphasized that the two-body bound states are characteristic to 
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FIG. 2: The two-electron scattering eigenfunctions for the incident two electrons coming in through 
the lead 1. The dotted circles surrounding the two electrons indicate the two-body bound states, 
which appear only when at least one of the electrons is reflected. 


open systems and do not appear under periodic boundary conditions. We also find that 
the two-body bound states are associated with the electron that is reflected at the QD. For 
example, if the incident two electrons come in through the lead 1 = £2 = 1 ), the two-body 

bound states appear only in the eigenfunctions g^il lk2 (x i,x 2 ), fi , i 2 1 fc 1 / C2 ( Xl >^ 2 ) and e fk 1 k 2 ( x ) 
since Z^ lk2 (x) = 0 as is depicted in Fig. [2j 

It is instructive to inspect the set of wave numbers characterizing each term of the eigen¬ 
functions in Eqs. (j22j) . As is found from the terms of the two-body bound states, the 
wave-number set {fci, /c 2 } of the incident states is not conserved and is scattered to the set 
{ed — W, E — ed + iT} including the imaginary part ffi. We note that the terms of the two- 
body bound states decay with the distance of the two electrons but are stationary in time 
since the total energy eigenvalue is real; the imaginary parts of the complex wave numbers 
cancel out each other in the total energy eigenvalue. By the completeness of the plane-wave 
functions |e 1 d Cl3;i+fc2:E2 )|fci, g M}, the terms are expanded as 


e Wf -ed)-r)bi-^2|+i£' £j 4 £2 



E-E' + iO 


c (A4,fc' 2 ;x 1 ,z 2 )e ,( *i*‘ + *i ra) , 


(25) 


where E' = k[ + k' 2 and c(k[, k 2 , x\, x 2 ) is the coefficient of the expansion. Hence we find 
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that the terms of the two-body bound states describe the various scattering processes to 
the sets { k \, k ' 2 } satisfying energy conservation E = E'. The coefficient c(k[, k' 2 ] X\, X 2 ) has 
the poles on the complex k\- and fc^-planes which come from the resonant pole k = £d — iT 
of the one-electron scattering eigenfunctions in Eqs. (fTTSlh Thus the two-body bound states 
appear as a consequence of one-body resonances^. 

The appearance of such many-body bound states is expected for general open QD sys¬ 
tems with localized interactions. We have shown that two-body bound states appear in the 


Anderson model with spin degrees of freedom^ and the double QD systems^. In Ref. [30, 


the A r -electron scattering matrix for another QD system with interactions was explicitly 
constructed in a real-time representation, where two-electron scattering eigenstates are ob¬ 
tained in an integral form. We speculate that, by evaluating the integral form, two-body 
bound states similar to ours should appear. 


C. iV-electron cases 

We can obtain the exact AT-clectron scattering eigenstates for arbitrary N. The scattering 
eigenstates for a few electrons can be constructed in a way similar to the two-electron cases. 
In the three-electron scattering eigenstates, for example, two electrons out of the three 
form the two-body bound states after the scattering at the QD= The explicit form of 
the scattering eigenstates for a few electrons leads to a conjectural form of the A’-electron 
scattering eigenstates. We have shown that they are indeed the eigenstates. 

We present only the results in the first order of Ug, which we need in the next section. 
The form of the A r -electron scattering eigenstates is given by 

\ E )= ('52 dNx 9rm-m N (Xl,-.., X N )J mi Oi) • • • C ] mN (: X N ) 

{mi} Jxi <- <XN 

+ 52 d N - 1 xe ni ... nN _ 1 (xi,...,x N - 1 )cl ll (x 1 )---cl lNl (x N - 1 )$)\0). (26) 

{m} 

Here we impose the antisymmetric relations for the A r -electron eigenfunctions as follows: 
9 m Ql -m QN (x Ql , ...,x Q n )= sgn (Q)g 

mp-mjv (oq,.. .,x N ), 

• -i x Rn- 1 ) = sgn(R)e ni ... riJV _ 1 (xi,... ,xjv-i), (27) 

where Q = (Q 1 ,... ,Qn) is a permutation of (1,2,..., N ) and R = (Ri ,..., Rn-i) is that of 
(1,2,..., A^ — 1). We consider the situation in which the electron with wave number /e, : , (?' = 
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1,... ,N) comes in through the lead li to the QD. The iV-electron scattering eigenfunctions 
are constructed in the first order of Ut as 




N 


E sgi 


i =1 


N—2 


Esgn(^) II 9 m P Q v k P ( x Qi) 9 ( x QN- 2 i ■ ■ ■ i X Q\) 


p,Q 


2 — 1 


Mpn-i^Pn') 


X tm Q N -i Uln Q N Z rnQN ,kp N _ lk p N 

0(U 2 ) : 


( x QN-iQN)e {kPN - 1+kPN)XQN °( x QNQN-i) e ( x QN-i ) 


(28) 


where P and Q are permutations of (1, 2,..., N), and 

( £i—£n) ( \ 


N -1 


Esgn(P) Yl 9n t P ,k Pi ( Xi ^ 


Mp N ) 


N J 
'N 


2 — 1 


TV-3 


i E sgn ( P ^) II ( x Qi) e fcpE 2 '^( x Q*- 3> • • •. x Qi) 


^PN- 2 ). 
‘n -2 


P,Q 


i=l 




N—2 


\ E sgn ( p< 5) n 


P,Q i=l 

tt y^ p N-P p N^ 

n QN-l Z ' n Q N _ 1 ^P N _ 1 k P N 


(“^Qiv-Je 1 ^- 1 +fcPjv )a:Qjv - 1 ^(^QAr-i) 


0(f/ 2 ), 


(29) 


where P is a permutation of (1, 2,..., IV) and Q is that of (1, 2,..., N — 1). Here we have 
used the notation 


d(x m , ...,x 2 , xi) = 6{x mim - 1 ) • • • 9(x 32 )0{x 2 i). (30) 

The third term in Eq. (1291) corresponds to the configuration in which one of the two electrons 
that form the two-body bound states is on the QD. On the other hand, the second term 
in Eq. (129|) corresponds to the configuration in which both of the two electrons that form 
the two-body bound states are in the leads, which has not appeared in the two-electron 
scattering eigenfunction ef^{x) in Eqs. (l22j) only in gfllk 2 i x i, x 2 ) hi Eqs. (l22]h In what 
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follows, we denote the eigenstate obtained by inserting the eigenfunctions in Eqs. (128)1 and 
(J2SD into Eq. ([26]) by \k\£) = |Aq ,... ... ,£ N ). 


IV. ELECTRIC CURRENT UNDER BIAS VOLTAGES 


A. IV-electron current 


By following the three steps of the extension of the Landauer formula given in Sec. Ill 61 
we next calculate the average electric current for the system under finite bias voltages^ 1 ^. 
First, we calculate the IV-electron current, that is, the quantum-mechanical expectation 
value of the electric-current operator / with the IV-electron scattering eigenstates | k\£). We 
assume ki ^ kj if £i = £j and restrict our calculation to the first order of Ue- We need to 
calculate the following overlap integral: 

(k]£\l\tj) = 2Im( / d N ~ l x ^ 

JxK -Kxn -i e{n . } 1 

(xi,... ,x JV _i,0)e ni ... njv _ 1 (xi,... ,xjv_i)), (31) 

where £ = ?> — £. By inserting the IV-electron eigenfunctions in Eqs. (T28|) and (129|) into 
Eq. fl3Tjh we obtain 


(k-,£\I\k-,£) 

\tit 2 \ 2 ( L \ N ~ l 


Tit 2 \2n / 
/ L \N-2 


£(-l) 4 Im(G fe ) 




\2ttJ 
2\t 1 t 2 \ 2 f L\N-a 


Tit 2 \2 ttJ 

0(U 2 




(32) 


Here L is the system length coming from the regularized square norm (k]£\k]£) = L/(2n) 
of the one-electron scattering eigenstates. In order to express the results, we have used the 
notation 


Gk 


1 

k - e d + ir ’ 


(33) 
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which is the one-electron Green’s function on the QD. We notice that the choice of the 
parameter a in Eqs. (J2J) simplifies the calculation. On the other hand, the square norm of 
the iV-electron eigenstates is calculated as 


(MIM) 


Y / d N x | g Ol, . . . ,Xat )| 2 + £ d X \c ni ...n N _i (xi, . . . , XjV— l) | 

fm Jxi< "<X N r„.l JX\<-<X N -! 


K}' 

= (s) a '- 2 (^)' v "£ £ ^(G*.*££«n-)«Sf^£K,(o))+o(^). (34) 


i^j m= 1,2 

It should be noted that the term in the (N — l)th order in L does not appear above. 
Combining Eqs. (I32j) and fjMjl . we obtain the N-e lectron current as 

(W\I\k;t) 


(k]£\k]£) 

_ jtiOp 2vr 
Tit 2 L 


N 


E(-l)'‘I ,n (GO 


2—1 


47T" 

T 2 


y fH w i«(0+)Z<S^(0) 


£(-l)'* £ U m 

m=l,2 

“TT'TT£ 


0{U 2 ). 




771—1,2 


(35) 


B. Average electric current 

Next, we take the reservoir limit of the iV-electron current in Eq. (l35|i to obtain the 
average electric current. We assume that the infinite lead substitutes for a large electron 
reservoir characterized by the Fermi distribution function / MjJ g(fc) = 1/(1 + e^ k ~^) with 
a chemical potential /i and an inverse temperature /?. We also assume that electrons are 
completely equilibrated in each electron reservoir before being re-emitted towards the QD, 
which is the main assumption of the extension of the Landauer formula. 

In the IV-electron current in Eq. (l35]h the A’-dependence appears only in the upper 
bounds of the sums on the wave numbers together with the factor 2tt/L, which means that 
we can take the reservoir limit Ng, L —» oo described in Sec. Ill E>[ It should be noted that 
the first term in Eq. (l35|i contains a single sum on i with the factor 2n jL and the second 
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term contains a double sum on i and j with (2tt/L) 2 while the third term is a double sum 
on i and j with (27r/L) 3 due to the square norm (k]£\k]£) appearing in the denominator. 
Therefore the third term in Eq. (ET)|) vanishes in the reservoir limit N#, L —y oo. 

In order to investigate the average electric current, we set i\ = • • • = = 1 and 

£n!+\ — • • ■ = — 2 in Eq. (155jl and relabel \k;£) by | k, h) with h z = k Nl+i , (1 < i < N 2 ). 

The ^-electron current is rewritten as follows: 


(k,h\I\k,h) 
(k, h\k , h) 


Jo + J\ + J 2 + 0(U 2 ), 


Jo = 
Ji = 


2^^ 27t 
7tT T 


JVi 


n 2 


E Im ( G d -Eim( G d 


i= 1 


i=1 


TiT 2 47t 
' tt 2 T L 2 


2 r 


u, 




v (t/iig;+t/rfg) +u 2 y. <£. 


( 22 ) 


l<i<N x 

1<J<JV 2 


l<i^j<N 2 


where we use 


(36a) 

(36b) 

(36c) 


A lm ) 

s kh 


Ill! 


G h 


1 )% m G k - (-1 ) m G h + 2i((-l ) l T m G* h G k - (-lr^GlG^ 


(37) 


and T i = |b| 2 /2. We omit the explicit form of J 2 in Eq. (I36a,p since it does not contribute 
to the average electric current in the reservoir limit L —y oo. Thus the parts of the 
iV-electron current that contribute to the average electric current are detemined by the 
two-electron scattering eigenstates. 

In the reservoir limit Ni, L —> oo, we replace the sums on ki in Eq. (136ap by the integral 
on k with / Mi/ g(fc) as 

2 N l roo 

— / dkf, lA (k)J(k), (38) 

i=1 7-A 

where we need to introduce the low-energy cutoff —A since the local dispersion relation of 
the lead is bottomless. At zero temperature (/^ = oo), the average electric current is given 
by 

(D = pHm( Gt ) - ^ E u, H k ( n h ^ + flh(2>) + 0(0*). (39) 

We notice that the first term in Eq. (j39l) reproduces the original Landauer formula in the 
noninteracting cases. The double summations in Eq. (I36cp give double integrals in the 
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second term in Eq. (1391) . which give a contribution of the Coulomb interactions. Through 
the integral formulas 


4 + 1 


/'/v 1 / -2 I 1 

I dkG k = ± log' * + 

/'/V 1 / 

/ d*G* = -(- ,, 

/—A TVe^-1 e A — 1/ 


+ i(arctan(e^) — arctan(e A )), 


1 

/ dkG* k Gk = — — (arctan(e^) — arctan(e A )), 
' —A r 


with = (ed — /x*)/T and e A = (ed + A)/T, we obtain the average electric current 


(40a) 

(40b) 


<4 = 


2^ 
7 TT 

rir 2 

tt 2 T 2 


-J- 


Uf ' ( F 4- “ r Ci) !°g 


1 = 1,2 


4 + ^ 

4 + i 


+ (T(q - e A ) - T*(2arctan(e*) - Ja))j 2 


0{U 2 


(41) 


where we use the notation 


j- = arctan(ei) — arctan(e 2 ), j\ = 2arctan(eA), 


,2 -s 


Js 


e 2 


4+i 4+i 


4 = 1 , 2 ). 


(42) 


We find that the average electric current (I) contains linear and logarithmic divergences in 
the limit A —y oo, which is similar to the symmetric case U\ = U^——. 


C. Universal electric current 


We employ a renormalization-group technique to deal with the divergences in the average 
electric current in Eq. (14111 . The divergences are due to the bottomless dispersion relation. 
By the renormalization-group analysis, we zoom into the Fermi energy and thereby discard 
all details that arise from the specifics of the dispersion relation. As a result, we obtain a 
universal form of the average electric current. 

We devise a Callan-Symanzik equatio n 24 ! 6 so that the average electric current may satisfy 
it. Let us introduce a parameter D = T^4 + 1. We can indeed see that, for D>T, |e d |, 
the average electric current (I) in Eq. (14TT) satisfies 


D 


_d_ 

dD 


+ 


£= 1,2 


d 

W 



(43) 
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where the beta functions /3r e and f3 ed are given in the first order of Ue as 

/Sr, = ~r e + 0(U 2 ), (3 ed = ——D + 0(U 2 ) (44) 

7T 7T 

with the average interaction U = (U d + U 2 )/2. The Callan-Symanzik equation of the 
form (1431) is an extension of the previous ones ——— for the case of the symmetric cou¬ 
plings. 

The general solution of the Callan-Symanzik equation determines a scaling form of the 
average electric current as 

(/) = D^T 2 , e d + , (45) 

where «/(-, ■) is an arbitrary three-variable function. Hence, if we change the parameters 

T e and e d as functions in D as 

T,0D) = T,(^) \ e d (D) = E d -^D (46) 

with the constants Tf, T 2 , (T K = T d + T 2 ) and E d , the average electric current (I) does not 
depend on D. The parameters T^(Z?) and e d (D) are referred to as renormalized parameters 
while the original parameters are called “bare” parameters, which we denote by T^ 0 , e dj o 
and D 0 . Here we fix the renormalized constants Ti and E d by the bare parameters as 

T t = E d = e d + -D 0 (47) 

\T k J 7r 

and express all physical quantities in terms of the renormalized ones. We shall see below that 
the sum T k = Ti + T 2 plays a role of a scaling parameter similar to the Kondo temperature. 

By inserting the renormalized parameters in Eqs. (ITSD into the average electric current 
(I) in Eq. (jH)) and rearranging it with respect to the interaction parameter Ui, we obtain 
the universal electric current 


(I) 


2T,T 2 , 
7tTk 


T{T 2 





+ 0 ( 7 / 2 ), 


Teji) log(e^ + 1) + (T K i ( - T t (2 arctan(e £ ) - tt ))j 2 

(48) 


where j_ and j s are j_ and j s with q = (E d — H()/Tk in place of q, respectively. Thus the 
average electric current (/), which was originally described by the bare parameters T^q, e d)0 , 
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Dq and Ue, is now characterized by the parameters Tg, E d and Ug. As a result, the linear and 
the logarithmic divergences of the average electric current are absorbed into the parameters 
Tn and E d . 

Let us consider the current-voltage (I-V) characteristics of the universal electric current. 
We put E d = 0 and consider the cases /p = — ji 2 — V/2 with the bias voltage V. Then we 
have 


T k 

4Ti r 2 


(/) 


1 

= — arctan 
7r 



1 

27T 2 


(U — 5) arctan 



+ 0(U 2 ), 


(■u + s ) 


VI (2T] 


k 


V 2 /(2T K y + lJ 


log 


V 


2Ti 


K 


+ 1 


(49) 


where 5 is the asymmetry parameter defined by 

«= gmLigmA . (so) 

27k ' ' 

We note that the parameters T x and T 2 depend on U\ and U 2 through Eqs. (J47]). 

We hnd from Eq. (1491) that the bias voltage V is scaled by the parameter Tk- In the case 
of symmetric interactions U\ = U 2 , we have 5 = 0 . Hence, after taking appropriate scaling 
factors for the electric current (I) and the bias voltage V, the I-V curve is independent of 
the parameters T { and T 2 . In the asymmetric cases with r 10 ^ r 2j o and U\ ^ U 2j on the 
other hand, the parameters T\ and T 2 play a nontrivial role since, even if we rescale the 
electric current, it still depends on 7\ and T 2 through the asymmetry parameter 5. 

We deal with the parameters T x and T 2 in the first order of U(. By solving the equation 
for Tk, which is obtained by the first equation in Eqs. (l47lh the parameter T \ is expanded 
in the hrst order of Ut as 


T e = W, 0 (l + ^log(^)) +0(U 2 ). (51) 

We remark that, for the consistency with the expansion, the ratio Dq/Yq should be restricted 
to the region Dq/Yq e 71 "/^. Then the asymmetry parameter 5 in Eq. (|50l) is expanded as 

s = {Ul ~ LAa ~ r ^) + o(^). ( 52 ) 

2T 0 

By a physical intuition, we expect 5 > 0 since the case U\ > U 2 should correspond to the 
case T 1)0 > T 2)0 . By expressing the parameters as r^o = T 0 (l ± q)/2 with 0 < 7 < 1 
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FIG. 3: I-V characteristics of the rescaled universal electric current for U = 0.5 with 5 = 
0,0.1,0.2,0.3,0.4. 


and U 1/2 = U(1 ± 7 ') with 0 < 7 ' < 1 in the cases with r 10 > r 2j0 > 0 and U\ > f/ 2 , 
we have S = U 77 ' + 0{U 2 ). Hence the asymmetry parameter <5 takes a value in the range 
0 < 6 < U + 0(U 2 ). The J-C curve of the universal electric current for U = 0.5 and 
0 < 5 < 0.5 is depicted in Fig. |3] 

We observe the suppression of the electric current for large bias voltages F > T K , that 
is, negative differential conductance. This is because the formation of the two-body bound 
states promotes the reflection of electrons at the QD, as is illustrated in Fig. [2j and the 
logarithmic term in Eq. ()49l) decreases the electric current. In the first order of Ue, the 
negative differential conductance shows the power-law behavior (I) oc (V/Tk)-^ u ~ s ' ) ^ where 
the asymmetry parameter 6 appears. This means that the suppressed electric current is 
restored by the asymmetry of the system parameters. In the case r 10 = r 2i o, we have 5 = 0 
in the first order of Ug and the universal electric current in Eq. (1491) depends only on the 
average interaction U. This is consistent with the Callan-Symanzik equation: by changing 
the variables and setting Ti = T 2 = T /2, the Callan-Symanzik equation (1431) is reduced to 


+ O= o 


d 


d 


(53) 


with the beta function fir = —UY/tt + 0{U 2 ). The general solution is given in the form 
(I) = J(D^Y,e^ + UD/n) with an arbitrary two-variable function J(-, ■). Then the negative 
differential conductance appears for U > 0, which is essentially the same as the previous 
results in the case of symmetric interactions U\ — 


Let us compare the present results with the renormalization-group (RG) results in Refs. 


31 
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and 


32 


in Ref. 
and 


the RG flow equations for the level width T were obtained in the second order of Uj 


32 


Although general local dispersion relations of the leads were adopted in Refs. |31 


32j . the details of the frequency dependence of the density of states of the leads did not 


play any role in their analysis. Indeed, the linear divergence in our average electric current 
(I) due to the linearized dispersion relations is removed by the RG technique with the 
Callan-Symanzik equation as we have described above. Their definition of the parameters 
Xi and T 2 , which was approximately derived from the RG flow equations for the level width 
T, is equivalent to ours in Eq. (14T1) in the first order of Up. 

We can confirm that the universal electric current in Eq. (1491) is consistent with that of 
the RG results^— in the first order of Up. By using the renormalized band width Yp(D) in 
Eq. (149j) . we have 


y 1 {d)y 2 {d) 


r(-D) 

TiT 2 r 

T k - 
1 r 
— 27T 


arctan 


V 


2r (D) 


arctan 


V 


2Ti 


K 


(U — 5 ) arctan 


V 


(U+5), 


V/(2T k ) 


lo g(^T 

1 K 


0(U 2 ). 


(54) 


2T k J v ~ ' ~'1/ 2 /(2Tk) 2 + 1- 

Then, by putting Yp = Yp(\V/2 + iTk|) and Y = Ti + T 2 , the universal electric current in 
Eq. (1491) is expressed by 

4f if 2 /V 


(I) = arctan ) + 0{U 2 ), 


v , , ,, (55) 

7rr V2T/ 

which is in the same form as the noninteracting cases. The expression in Eq. (1551) agrees 
with that obtained in the RG results (see Appendix [A|. 
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We remark that, although the expression of the universal electric current in Refs. 


31 


and 


agrees with ours in Eq. (|55l) . their treatment of the renormalized parameters T\ and T 2 
included in the universal electric current is different from ours; they used the higher-order 
terms of Up in the defining relations of and T 2 in Eq. (I47]i while we have treated them in 
the first order of Up as is given Eq. (15TT) . As a result, even in the case T 10 = r 2 0 of symmetric 
lead-dot couplings, they observed the restoration of the suppressed electric current. We have 
shown that, for the asymmetry parameter <5 treated in the first order of Up, the restoration 
due to the asymmetric interactions does not occur at T^o = T 2 i o. In other words, in this 
case, there should be no restoration of the suppressed electric current for small Up, which 


seems to differ from the results of Refs. 


31 


and 


32 
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V. CONCLUDING REMARKS 


We have studied the average electric current for the open QD systems described by the 
two-lead IRLM in the asymmetric settings. By using the extension of the Landauer for¬ 
mula with the many-electron scattering eigenstates, we have calculated the average electric 
current for the systems under finite bias voltages. The calculation is in the first order of 
the interaction parameters, but otherwise we have not employed any approximations. The 
calculation itself has been considerably simplified compared to the previous one^^ treating 
the symmetric cases with the even-odd transformation. 

Through the renormalization-group technique with the Callan-Symanzik equation, we 
have obtained the universal electric current characterized by the scaling parameter TR and 
the asymmetry parameter S. The Coulomb interactions around the QD give rise to the 
negative differential conductance through the formation of the two-body bound states, which 
is a new point of view clarified in our analysis. Through the investigation of the asymmetry 
parameter <5, we have confirmed in the first order of U\ and U 2 that the suppressed electric 
current is restored in the asymmetric cases satisfying both T 10 7 ^ T 20 and U\ ^ U 2 . 

The analytic form of the universal electric current has enabled us to compare it with 
those of other approaches correctly. Our universal electric current has the same functional 
form as that obtained with the RG approach^^; it also reproduces previous result a 2 ° i23 ~ — 
in the symmetric cases. However, for U\ 7 ^ U 2 and T 1)0 = T 2 j o, our results indicate that there 


is no restoration of the suppressed electric current to first order in U\ and U 2 , while Refs. 


31 


and 


32 


indicate that there is. This restoration may result from higher-order terms in U\ and 


U 2 . To verify its validity analytically, we need a consistent treatment of these higher-order 
terms, which can be done with our extension of the Landauer formula. 

The key element for the practical calculations in the extension of the Landauer formula 
is the explicit form of many-electron scattering eigenstates. The present calculation in the 
first order in the interaction parameters can be extended to higher orders by using the 
exact 7V-electron scattering eigenstates that we have already obtained. We expect that such 
calculation should be applied to other physical quantities of the open QD systems such as 
the dot occupancy. 
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Appendix A: Results of the RG flow equation for T 


In our settings, the universal electric current obtained in Refs. [3l| and [32] is expressed as 

2f 1 f 2 


(-Org = 


7lT 


^i-e d \ 7^2-ed 

arctan —-— — arctan 


(Al) 


r / v r 

where is the solution of the RG flow equations for the level width of the QD and T = 
Ti + T 2 . It was approximately determined by the self-consistent equation 

D 0 


r f ~r 


1,0 
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|T - i(n e - e d )| 

with gi = Uij7 T + 0{U 2 ). In the first order of U(_, the approximate solution is given by 

T k 


(A2) 


Y f — Tf 


= Tt 


|f - i (n e - e d )| 


u± 


+ 0(U 2 ) 


- K 


u t 


vr < yy + °(U 2 
\T k - - e d )| / 

= r,(|(/i,-e d ) + iT K |) + 0(f/ 2 ), 


(A3) 


where Y^D) is the renormalized level width defined in Eq. (I46lh By inserting this into 
Eq. (lAlj) and setting /ii/ 2 = ±.V/2 and e d = 0, the expression in Eq. (IA1D agrees with 
Eq. (J55D. 
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